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ABSTRACT: In molecular vibrational infrared spectroscopy, absorption spectra arise from molecular 
vibration and correspond to transitions between the vibrational energy levels associated with a given 
electronic state of the molecule. The vibrational transitions, which fall in the near infrared region, are 
induced through the interaction of the molecular electric dipole with the electric vector of the 
electromagnetic radiation. The near infrared region extends roughly from 1umto 10?um. The article 
explains the pure vibrational absorption spectra of diatomic molecules such as HCl, HBr, HI, CO, ... etc. In 
order to explain the vibrational spectra, diatomic molecules are treated as harmonic oscillator and 
anharmonic oscillator. In the harmonic oscillator model, we get only one absorption band at the 
wavenumber value@,,, corresponding to frequency of oscillationv,,, while in the actual experimental data, 
there are many absorption bands corresponding to wave numbers slightly lesser than @ 5.) 2Wosc, 3Wosc» 
....1 he occurrence of these additional bandsis attributed to the selection rule Av = +2, +3, +4, ...... The 
additional bands are having lesser intensity and are called overtone bands. 


KEYWORDS: Pure vibrational spectra, Diatomic molecules, Harmonic oscillator, anharmonic oscillator, 


Population of excited states, Fundamental band, Overtone bands, Hot band. 


https: //doi.org/10.29294 /IJASE.8.1.2021.2016-2024 


1.INTRODUCTION 

Infrared (IR) spectroscopy, an analytical 
technique that takes advantage of the vibrational 
transitions of a molecule, has been of great 
significance to scientific researchers in many fields 
such as_ protein characterization, nanoscale 
semiconductor analysis and space exploration [1]. 
Infrared (IR) spectroscopy is one of the most 
common and_ widely’ used _ spectroscopic 
techniques employed mainly by inorganic and 
organic chemists due to its usefulness in 
determining structures of compounds and 
identifying them [2]. Chemical compounds have 
different chemical properties due to the presence 
of different functional groups [3]. 

The atoms in a molecule do not remain in fixed 
relative position but vibrate about some mean 
position. The internuclear distance changes when 
the nuclei of these molecules vibrate relative to 
each other. Thus, these molecules possess an 
oscillating dipole, which according to classical 
electrodynamics, emits radiations of frequency 
that fall in near infrared region [4]. Heteronuclear 
diatomic molecules such as HCl, HI, CO, CN, .... etc. 
have intrinsic electric dipole moment and are 


capable of interacting with electromagnetic 
radiation. Hence, these molecules’ exhibit 
vibrational spectra. Homonuclear diatomic 


molecules do not have dipole moment and hence 
do not give vibrational spectra [5, 6]. 
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To develop a description of the vibrational states, 
we will consider the diatomic molecule to be a 
Harmonic oscillator, i.e. atoms execute simple 
harmonic motion about their mean position and 
the restoring force applied on the atoms is directly 
proportional to displacement. A _ theoretical 
analysis of the quantized vibrational energy levels 
and the vibrational transitions exhibited by the 
polar heteronuclear molecules is presented. The 
harmonic oscillator model is a_ great 
approximation of a molecular vibration, but has 
key limitations:1. Due to equal spacing of energy, 
all transitions occur at the same frequency (i.e. 
single line spectrum). However, experimentally 
many other lines (called overtones) are often 
observed 2. The harmonic oscillator does not 
predict bond dissociation; you cannot break it no 
matter how much energy is introduced. The 
discrepancies of theoretical and experimental data 
have been resolved by considering the diatomic 
molecule as an anharmonic oscillator. 


2. THEORY OF PURE VIBRATIONAL SPECTRA 


Now, we will discuss the theory of vibrational 
spectra of diatomic molecules. There are two 
models to discuss the pure vibrational motion of 
heteronuclear diatomic molecule. 
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i. Harmonic oscillator model 
ii. Anharmonic oscillator model 


2.1 DIATOMIC MOLECULE AS A HARMONIC 
OSCILLATOR 


Consider a diatomic molecule consisting of two 
atoms Aand B and separated by equilibrium 
distance r,, .Fig. 1 shows how the potential energy 
of diatomic molecule varies with the internuclear 
distance r. The minimum of this curve 


Energy 
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corresponds to the normal configuration of the 
molecule. The zero of the curves is found at the 
value r=r,,. Any energy in excess of this, for 
example, ¢€, arises because of extension or 
compression of the bond [1]. If one of the atoms 
(suppose A) is considered to be stationary on the 
r = 0 axis, the other will oscillate between B’ and 
B” If the energy is increased to €,, the oscillation 
will become more vigorous [1]. 


internuclear distance 


Fig 1. Parabolic curve of energy plotted against the expansion and compression of harmonic 
oscillator obeying Hooke’s law 


Suppose the bond is distorted from _ its 
equilibrium length 7, to new length r. The bond, 
like spring, obeys Hooke’s law. Then 


f= -k(r- 1%.) -() 


Here, f is the restoring force, k is the force 
constant, and r the interaction distance [3]. 

In this case, the potential energy curve is 
parabolic and has the form 


1 2 
V= 5kr - Teg) . (2) 


This model of vibrating diatomic molecule is 
called “Simple Harmonic Oscillator” model [1, 2]. 
Fig. 1 shows the parabolic curve of energy plotted 
against the expansion and compression of 
harmonic oscillator. 

Classically it is easy to show that the oscillation 
frequency is 


1 jk 
Vosc == |{—--- (3) 
27 | 


where, is the reduced mass of the diatomic 
molecule. 

The energy of Harmonic oscillator is given by 
[after solving Schrodinger equation] 


1 
Ey = Avosc (v + 5) (4) 


where, v is called the vibrational quantum number 
which can take values: v = 0, 1, 2, 3,........ Equation 
(4) gives the allowed energies for the harmonic 
oscillator. 

A special feature of the quantum mechanical 
oscillator is the existence of zero-point energy, 


1. 
7 AVosc: 
Let us now investigate the expected spectrum of 


such an oscillator. The vibrational terms (energies 
in wavenumber units m1 or cm~') are 
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E, Vosc 
G — 
aaa 


(v +3) (5) 


Vv 4 7 . 
where, “~ is the classical frequency in wave 


number unit and is called as vibrational constant 
denoted by W,5c. 
Thus, 


G(V) = Dose (v + 5) ... (6) 


Substituting v = 0, 1, 2, 3,......... , we get 
1_ 
G(0) = 7 ose 
3 
G(1) = 2 Pose 
5 
G(2) = 2 Vosc 


= 
G(3) = 2 Vosc 


v=6 
v=5 
Ss v=4 
a | 
EF v=3 
= 
— 
v=2 
v=1 
v=0 
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Thus, we have a series of equidistant discrete 
vibrational levels, the common spacing being 
@osc[1]. 

When a transition takes place between a lower 
level v' and an upper level v’, the wave number of 
the absorbed electromagnetic radiation is 


Veo! = G(v’) _ Gv) 


= _ nv 1 _ oe | 
Vy’ sy" = ®osc (v + 5) — Wosc (v + >) (7) 
Now, for a harmonically oscillating molecule of 
dipole momentM, the matrix element is found to 


be non-vanishing only when 


on 0 andA +1 (8) 
—+#0andAv=+ ive 
dx 
1 
I 
1 
i 
13 — 
2 Cae 
Als... 
2 
a 
zs Pose. 
= amie: 
yA 
5 
= Dose. 
2 
a 
2 
1A Dose 
2 


= __. 


DWosc. 


nr 


Wave number V (cr) 


Fig 2. The allowed vibrational energy levels and transitions between them for a diatomic molecule 
undergoing simple harmonic oscillation 


That is, vibrational transition can only occur 
when the molecule has a permanent dipole 
moment M which changes’ with distance x. 
Further, the selection rule Av = +1 gives [1] 

vi=vt¢1 (-v'>Vv)..Q9) 

The spacing between the vibrational energy 
levels is considerably larger than the spacing 


between the rotational energy levels of a molecule, 
in fact larger than k,T at room temperature. 
Hence, most of the molecules in a sample exists in 
the v = 0 state, with only their zero-point energies 
(except in case of very heavy molecules which 
have smaller vibrational spacings @,,,). 

Thus, the main vibrational transition in 
absorption is vy =0>v'=1. 
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The absorbed wave number corresponding to the 
transition (v’ = 0tov" = 1) is 


2 = Get. ee ads 
Vy' ov" = osc (v + 5) — ®osc (v + 5) 


a = 1 = 1 
Vos1 = Wosc (1 + 5) — Wosc (0 + 5) 
_ = ae 

Vos1 = 7 os _ 7 ose 


c 
Vo1 = Dosc + (10) 


Thus, the vibrational spectrum is expected to 
consist of a single absorption band at @,,,. Fig. 2 
shows the energy level diagram with transitions in 
a diatomic molecule undergoing harmonic 
oscillations. The observation of a very intense 
band in the infrared spectrum thus leads to the 
conclusion that it is a vibrational spectrum, the 
nuclei carrying out harmonic vibrations along the 
internuclear axis [1]. 


However, the occurrence of weak overtone 
bands at wave numbers nearly two, three times 
the wave numbers of the main band is still to be 
explained, and is connected with the deviations 
from the oscillations being “harmonic” [1]. 


2.2 THE MOLECULE AS ANHARMONIC 
OSCILLATOR 


A comparison of an observed near infrared 
spectrum with that expected from a diatomic 
molecule treated as harmonic oscillator shows an 
important disagreement. The harmonic oscillator 
would give a single band at wave number @,,,, 
which is the classical frequency of vibration of the 
molecule. The actual infrared spectrum is, 
however, found to consist of an_ intense 
(fundamental) band at @,,,, plus a number of 
weak bands (overtones) at wave numbers slightly 
lesser and lesser than 20 5c, 3@osc) vse The 
observation indicates that the selection rule Av = 
+1 is not strictly obeyed [2, 3], and transitions 
corresponding to Av > 1 do take place. This, in 
turn, is attributed to the fact that the dipole 
moment of the molecule is not strictly linear with 
respect to the internuclear displacement x (= r — 
Toq.): This is expressed as _ “electrical 
anharmonicity” of the molecule [7-9]. The 
observation that the overtones appear not exactly 
at 2Wosc1 3Woscy vere but at lesser and lesser 
values indicates that the vibrational energy levels 
are not exactly equally spaced, but converge 
slowly. This is attributed to the fact that for an 
actual molecule the potential energy curve is not 
strictly parabolic (except near the minimum). That 
is, the potential energy function V(r) is not 
harmonic and we must include terms higher than 
quadratic in the Taylor’s series expansion of V(r). 
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This is expressed as “mechanical anharmonicity” 
of the molecule [2, 8]. 

A purely empirical expression which fits the 
potential energy curve (observed) to a good 
approximation was derived by P.M. Morse, and is 
called the Morse equation [1, 7]: 

V(r) = D, [1 — ef ar-rea)}]” ... (11) 
where, a is a constant for particular molecule and 
D, is the depth of the potential energy curve. 
To a first approximation, let us include the square 
and cubic terms in V(r), so that 


V(r) = f(r eq.) — g(r ~ Teg) (42) 
where g < f. 


Substituting the value of V(r)from equation (13) 
in the Schrodinger equation and solving by 
perturbation method, we get the energy eigen 
values of the anharmonic oscillator as 


_ ih 2 1)? 
E,, = hc@, (v + 5) —hceWeX¢e (v + 5) 
+ ikante (13) 


The corresponding term values are given by 


ae 
ve Be ; 


G(v) = @. (v + 5) — DeXe (v + 5) (14) 


The quantity @, is the wavenumber spacing of 
energy levels that would occur if potential energy 
curve were a parabola, @,x, is the “anharmonicity 
constant” which is much smaller than @, and is 
always positive [4, 6]. 

The equation (14) shows that the energy levels 
of anharmonic oscillator are not equidistant, but 
their separation decreases slowly’ with 
increasing v. Fig. 3 shows the energy level diagram 
of a diatomic molecule undergoing anharmonic 
oscillations. 


Substituting v = 0,1, 2, 3,...00 , we get 

1 1 

G(0) > 7 We a 4 WeXe 
3. 9 

G(1) = 57 We — G Vere 
5. O25 

G(2) = 2 le 4 Vere 
7 49 _ 

G(3) = 2 le a ee 


When a transition takes place between a lower 
level v' and an upper level v’ (due to absorption), 
the wave number of the emitted radiation is 


Vga" = G(v’) ~ Gv) 
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Vy’ sy" = @(v" - v') -— &,x.(v" -v’)(v" + v' + 1) 


vo 


Thus, as v increases, the separation between 
successive bands (or levels) decreases very nearly 
linearly, in agreement with observation [2, 5]. 


Let us now investigate the infrared spectrum for 
the anharmonic oscillator. Since the eigen function 


Energy 
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of anharmonic oscillator are very similar to those 
of the harmonic oscillator, the selection rule Av = 
+1 still holds giving the most intense transitions. 
In addition, for the anharmonic oscillator, 
transitions corresponding to Av=+2, +3, 
ued also appear, even though with rapidly 
decreasing intensity. The possible transitions in 
absorption when all the molecules are initially in 
the v = 0 state are shown in Fig. 4. This explains 
the appearance of observed weak overtone bands 
together with intense fundamental band [1, 10]. 

It is also seen that transitions with Av = 2, 3, 
4,..... have approximately, but not exactly, two, 
three, four, ....... times the wave number of the 
transition Av = 1, in agreement with observation. 


internuclear distance 


Fig 3. The Morse curve: the vibrational energy levels of diatomic molecule undergoing anharmonic 
expansions and compressions 


1. Fundamental Band: The transition from v’ = 0 
to v’ =1 gives most intense fundamental band 
with wavenumber [4, 5] 

Vo +1 = (1 — 0) — @, x, (1 — 0)(1 +0 + 1) 


Vo +1 = We — 2@,X%- = W(1 — 2x,)... (16) 


2. First Overtone Band: The transition from v’ = 

0 to v = 2 gives weak band (with lesser intensity 

than fundamental band) with wavenumber [1, 5] 
Vp +2 = (2 — 0) — B x,(2 — 0)(2 +0 + 1) 


Vo +2 = 20, — 60,X, = 2@,(1 — 3x,) ... (17) 
3. Second Overtone Band: The transition from 
v =0 to v =3 gives much weaker band with 
wavenumber 


Ip 53 = B(3 — 0) — &x-(3 — 0)(3 +0 + 1) 
Ip +3 = 3W, — 120,x, = 3@,(1 — 4x,) ... (18) 


4. Third Overtone Band: 
The transition from v =0 to v =4 gives 
muchweaker band (negligible intensity) with 
wavenumber 

Vo +4 = We(4 — 0) — @,.x.(4-0)44+0+4 1) 


Ip +4 = 40, — 200,x, = 40, (1 — 5x,) ... (19) 


Thus, the overtone bands are observed at wave 
numbers slightly lesser and lesser than 2@,, 3@,, 
AD veereee (See Fig. 4). 
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Energy 


E-ISSN: 2349 5359; P-ISSN: 2454-9967 


Fig 4. The allowed vibrational energy levels and some transitions between them for a diatomic 
molecule undergoing anharmonic oscillations 


3. ISOTOPE EFFECT ON VIBRATIONAL ENERGY 
LEVELS 


Different isotopic molecules have different 
vibrational levels, and hence different vibrational 
frequencies. From equation (3), the classical 
frequency of a molecule (assumed as harmonic 
oscillator) is given by [1] 


_ 1 (fk 
~~ 2ae |e 


Vosc 


The force constant k is determined by the 
electronic motion only and is therefore exactly the 
same for different isotopic molecules. The reduced 
mass is, however, different for different isotopes. 
If w' is vibrational constant for the heavier 
isotope, we have [1] 


ow! ee u 
w Vosc u 


4. THERMAL DISTRIBUTION OF VIBRATIONAL 
ENERGY LEVELS (POPULATION RATIO) 


According to Maxwell-Boltzmann distribution 
law [1], the number of molecules in the vth state, 
N,, relative to that in the lowest state, No is given 


by 


N= Nye~ Ev—F0)/keT 


where, E,, = hc@, (v + 4) — hcWeXe (v + 5) = 
hceWe (v + 5) 

(ignoring the second term with anharmonicity) 

Ny 


— = e~(Ev-Eo)/kaT | (21) 
) 


Now, consider the spacing between the 
vibrational energy level as 1000 cm, ie.@, ~ 
®osc = 1000 cm~?. Let us consider the population 
of ground state is Np. 


Now, we will calculate and compare the 
population of first excited state, second excited 
state, third excited state relative to ground state. 


(i) Ratio of populations of first excited state 
and ground state: 


N. _C1-Eo) - (*2ase) 
l=, kT =e \ keT 


N, 6.626 x 10-34 x 3 x 102° x 1000 
" 1.38 x 10-23 x 300 


= exp(—4.78) = 0.00822 
~ 0.01 0r1% 


N, = 1% ofNy 


(ii) Ratio of populations of second excited state 
and ground state: 
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(E2-Eo) (7eeese) 
Nz keT =e \ kpT 


Np _ 2 x 6.626 x 10734 x 3 x 10° x 1000 
i 1.38 x 10-23 x 300 


= 0.0000675 ~ 0.00007 or 0.007 % 
N, © 0.007 % ofN, 


(iii) Ratio of populations of third excited state 
and ground state: 


N. (E3-Eo) (MeBesc) 
3 — 9 keBT =e \ kpT 


3 x 6.626 x 10734 x 3 x 107° x 1000 
1.38 x 10-23 x 300 
= 0.000000555 
= 0.0000006 or 0.00006 % 


Nz © 0.00006 % of No 


From above calculations, it is clear that the 
population of the first excited state(v = 1)is 
approximately 1% of the population of ground 
state(v = 0), the population of the second excited 
state (v = 2) is approximately 0.007 % of the 
population of ground state(v=0) and the 
population of third excited state (v=3) is 
0.00006% of the population of ground 
state(v = 0). Thus, at low temperature we may 
ignore all transitions originating from vibrational 
quantum numbers v > 1 and restrict ourselves to 
the transitions originating at v = 0 only [1]. 


Effect of increasing temperature from T = 
300K to 600K (at @,,, = 1000cm~1) 


Ratio of population of molecules in the first 
excited state to the ground state is 


N. (E1-Eo) (“ese 
ah ag ket =e \ aT 


1.38 x 10-23 x 600 
= exp(—2.40072464) = 0.0907 


N, ( 6.626 x 10734 x 3 x 107° x ) 
— xX CC ooo 


N, ® 9.07 % ofNo 


Here, the population of the first excited 
state(v = 1) is approximately9 %of the 
population of ground state(v = 0). Thus, at T = 


5. DISSOCIATION ENERGY OF DIATOMIC 
MOLECULE 

The energy of vibrational energy level may be 
written as: 
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600 K, we may consider transitions originating 
fromv = 1. The transitions starting from v = 1 are 
known as hot bands. 


Effect of decreasing vibrational frequency from 
1000cm~1 to 500cm“ (T = 300K): 


(E1-Eo) he®osc 
a ee te 
No 
Ny _ 6.6 x 10734 x 3 x 107° x 500 
No P 1.38 x 10-23 x 300 


= exp(—2.40072464) = 0.0907 
N, ® 9.07 % of No 


In other words, the population of the v = 1 state 
is nearly 9 percent of the ground state population. 
Thus, we can say that a similar increase in the 
excited state population would arise if the 
vibrational frequency is lowered from 1000 cm 
to 500 cm7?. Thus, at T = 300 K, we may consider 
transitions originating at v = 1 if the vibrational 
frequency is lowered. 


Hot Bands: The transitions from first excited state 
(v=1) to the higher states (v>1) are not 
observed unless the temperature of the molecular 
sample is high or if the vibration has a particularly 
low frequency. This is due to the fact that the 
population in v = 1state is negligible at lower 
temperature and it becomes appreciable only at 
higher temperatures [1]. 


The transition v’= 1 — v” = 2 appears at wave 
number 

Vy’ py" = De(v" — v’) — O,x.(v" — v)(v" + v' +1) 
Trop = Dp(2 — 1) — Ox,(2-D(2+1+1) 

Tyo2 = Do — 4D oXe 


452 = @(1 — 4x,) ... (22) 


The transition v’=1— v”=2 is very weak 
and it will be observed closer to and at slightly 
lower wave numbers than the fundamental band 
[1]. These transitions are generallyknown as ‘hot 
bands’, since high temperature is the one 
condition for their existence. Their nature may be 
confirmed by raising the temperature of the 
sample when a true hot band will increase in 
intensity. 

2 


1 1 
E,, = hc@, (v + 5) —hcWeXe (v + 5) .. (23) 


And therefore, separation between neighbouring 
levels, AE, is normally 


AE = Fyat ~ Ey 
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= 1 = 1\2 
AE= [aca (v+1++) —ncw,x, (v+1++) |- 
= 1 = 1\2 
[aca (v + *) — hce@eX¢e (v + *) | 
AE = hc@, — 2hc@,.x,(v + 1) 


AE = hco,[1 — 2x,.(v + 1)] 


This separation obviously decreases with 
increasing v and the dissociation limit is reached 


2023 
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it 7 1 ty 
Ema = NEW, ga) 


Eas = hee gE] ~ h(E) 
max = Me Va al \ax, 


This will represent the dissociation energy of the 
molecule with respect to the bottom of the 
potential energy diagram. The true dissociation 
energy will be given by 


when AE — 0. Thus, the maximum value of v is Do = De — Ey ao = heii, ae 
given byv,,,,,, where F 4x, ; 4 
AE = hea, [1 — 2%2(Vmax + 1)] = 0 ine [5a = jae%| 


[1 — 2x.(Umax + 1)] = 0 


meh Xe 1 Xe 
Eee ly he 
Vmax = 2x. -1 ies (24) 1 
e Do = hc@, 4x, | as (25) 
If we recall that the anharmonicity constant, 
For HCl molecule, @, =2886cm 1 and 


x. ~1073; hence vygy = —— —1=500-1= dissociation energy (for x, = 0.001) 

499 =~ 500. 1 1 
ee 

Substituting value of v,,,, in equation (1), we get 4x, 2 


1 1)? 
Emax = De = hcWe (Gt + 5) — ACWeX-e Co + 5) 


= (6.626 x 10-34) x (3 x 108) 
x (2886 x 100) 


mi 2) 
x ———— se 

4x0.001 2 

= 1431329304.6 x 10-7°Joule 

= 0.89458 eV = 0.9 eV 


Thus, the dissociation energy of HCl molecule is 
approximately 0.9 eV. 


TABLE 1: The population ratio of first, second, third, fourth and fifth excited state with the ground 
state at 300 K, 600 K and 900 K. 


Energy N,/No N,/No N,/No 
Level /state (atT = 300K) | (atT = 600K) (at T = 900K) 
First Excited State 0.00822 0.0907 0.202 
Second Excited State 0.0000675 0.00822 0.0407 
Third Excited State 0.000000555 0.000745 0.00822 
Fourth Excited State 0.00000000456 0.0000675 0.00166 
Fifth Excited State 0.00000000003' 0.00000612 0.000335 
Rohit Singh 


www.mahendrapublications.com 


International Journal of Advanced Science and Engineering 


Int. J. Adv. Sci. Eng. Vol.8 No.1 2016-2024 (2021) 


6. CONCLUSIONS 


This article explains the infrared spectra of 
heteronuclear diatomic molecules such as HCl, 
HBr, HI, CO, ... etc. In order to explain the 
vibrational spectra, diatomic molecules are 
treated as harmonic oscillator and anharmonic 
oscillator. In the harmonic oscillator model 
infrared spectra are very simple; only the 
fundamental transitions, Av = +1, are allowed 
while in the actual data, there wereadditional 
bandspresent corresponding to wave numberclose 
tO 2Wosc. 3Woscr 4Woscs The additional bands 
are having lesser intensity and are called as 
overtone bands. The occurrence of these 
additional bands is attributed to the selection rule 
Av = +2, +3, +4, In this article, we have also 
explained the effect of isotopic substitution on the 
infrared spectroscopy. The expression of 
dissociation energy of molecules is also discussed. 
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